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Abstract 

Limits on the cross-section times branching fraction for neutral Higgs bosons, pro- 
duced in pp collisions at ^/s = 7 TeV, and decaying to two tau leptons with pseu- 
dorapidities between 2.0 and 4.5, are presented. The result is based on a dataset, 
corresponding to an integrated luminosity of 1.0 fb~^, collected with the LHCb de- 
tector. Candidates are identified by reconstructing final states with two muons, a 
muon and an electron, a muon and a hadron, or an electron and a hadron. A model 
independent upper limit at the 95% confidence level is set on a neutral Higgs boson 
cross-section times branching fraction. It varies from 8.6 pb for a Higgs boson mass 
of 90 GeV to 0.7 pb for a Higgs boson mass of 250 GeV, and is compared to the 
Standard Model expectation. An upper limit on tan /3 in the Minimal Supersym- 
metric Model is set in the m'^o'^ scenario. It ranges from 34 for a CP-odd Higgs 
boson mass of 90 GeV to 70 for a pseudo-scalar Higgs boson mass of 140 GeV. 
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1 Introduction 



The discovery of a boson with a mass of about 125 GeV by the ATLAS [1] and CMS [2] 
collaborations requires further investigations to confirm whether its properties are com- 
patible with a Standard Model (SM) Higgs boson or if it is better described by theories 
beyond the SM, such as supersymmetry. The ATLAS and CMS measurements have been 
made at central values of pseudorapidity, rj; investigations in the forward region can be 
provided by the LHCb experiment, which is fully instrumented between 2 < 77 < 5. Both 
measurements of cross-sections and branching fractions allow different models to be tested. 
In this paper, model- independent limits on the Higgs boson-f cross-section times branching 
fraction into two tau Icptons arc presented for the forward region and compared to SM 
Higgs boson predictions. Model-dependent limits for the Minimal Supersymmetric Model 
(MSSM) Higgs bosons, in the scenario where the lightest supersymmetric Higgs boson 
mass is maximal (m™o^^) [3], are also given for the ratio between up- and down- type Higgs 
vacuum expectation values (tan /3) as a function of the CP-odd Higgs boson {A^) mass. 

2 Detector and datasets 

The LHCb detector [4] is a single-arm forward spectrometer. The components of par- 
ticular relevance for this analysis are a high-precision tracking system consisting of a 
silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon- 
strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, 
and three stations of silicon-strip detectors and straw drift tubes placed downstream of the 
magnet. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. The trigger [5] consists of a hard- 
ware stage, based on information from the calorimeter and muon systems, followed by a 
software stage, which applies a full event reconstruction. 

Simulated data samples are used to calculate signal and background contributions, 
determine efficiencies, and estimate systematic uncertainties. Each sample was generated 
as described in Ref. [6], with Pythia 6.4 [7] using the CTEQ6L1 leading-order PDF 
set [8] and passed through a Geant4 [9] based simulation of the detector [10]. The 
LHCb reconstruction software [11] was used to perform trigger emulation and full event 
reconstruction. 

The dataset used for this analysis is identical to that described in our previous mea- 
surement of the Z cross-section using tau final states [12], which corresponded to an 
integrated luminosity of 1028 ± 36 pb~^, taken at a centre-of-mass energy of 7 TeV. The 
Z TT decays are identified in five categories: t^t^, T^Tg, TgT^, t^t/j and T^Th, defined so 
as to be exclusive, where the subscripts indicate tau decays containing a muon (//), elec- 

*The symbol $° is used throughout to indicate any neutral Higgs boson. Additionally, charge conju- 
gation is implied and the speed of light is taken as 1. 
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tron (e), or hadron (/i) and the ordering specifies the first and second tau decay product 
on which different requirements are applied. The first tau decay product is required to 
have transverse momentum, px, above 20 GeV and the second to have px > 5 GeV. Both 
tracks are required to have pseudorapidities between 2.0 and 4.5, to be isolated with little 
surrounding activity, to be approximately back-to-back in the azimuthal coordinate, and 
their combined invariant mass must be greater than 20 GeV. The tracks in the t^t^, T^Th, 
and TeT^ categories are required to be displaced from the primary vertex. Additionally, 
the TfjTn category requires a difference between the px of the two tracks and excludes 
di-muon invariant masses between 80 and 100 GeV, to suppress the direct decays of Z 
bosons into two muons. Full details on the selection criteria can be found in Ref. [12]. 

The invariant mass distribution of the two final state particles for the selected $° tt 
candidates is plotted in Fig. 1 for each of the five categories separately and combined 
together. No candidates are observed with a mass above 120 GeV. The distributions of 
Fig. 1 differ from those of Ref. [12] as the simulated mass shapes arc calibrated to correct 
for differences between data and simulation, and the Z ^ tt distributions are normalised 
to theory. 

Six background components are considered: Z — > tt; hadronic processes (QCD); elec- 
troweak (EWK), where one r decay product candidate originates from a W or Z boson 
and the other comes from the underlying event; ti; WW; and Z ^ i£ where ££ indicates 
electrons or muons originating from a leptonic Z decay. 

All backgrounds, except Z — )• tt, have been estimated in Ref. [12]. The distribution 
and normalisation of QCD background events is found from data using same-sign events. 
The electroweak invariant mass distribution is taken from simulation and normalised using 
data. The small contributions from ti and WW production are taken from simulation, 
while the Z ii invariant mass shape and normalisation are determined from data. 

The invariant mass distributions for $° tt and Z ^ tt decays are evaluated from 
simulation where the mass resolution has been calibrated using the Z ^ mi invariant 
mass peak. Each event is re- weighted by a factor (u x £)/(crsiin x Sgij^), which provides 
a negligible correction in comparison to the mass rcsohition calibration. The efficiency, 
e, for triggering, reconstructing and selecting candidates has been evaluated as a func- 
tion of momentum and pseudorapidity using data-driven techniques and is described in 
Ref. [12], while e^^^^ is the corresponding efficiency in simulation. The cross-section for the 
process in simulation is represented by (Tsim, while a is the theoretical cross-section. The 
Z ^ TT sample is normalised using the cross-section calculated with Dynnlo [13] using 
the MSTW08 PDF set [14]. The tt signal distribution is found from simulated 

gluon-fusion events. The signal samples were generated in mass steps of 10 GeV from 
90 GeV to 250 GeV. For both the SM and MSSM Higgs bosons, the normalisation of the 
signal uses the theoretical calculations described below. 

The SM cross-sections, using the recommendations of Refs. [15] and [16], are calcu- 
lated at ^/s = 7 TeV with the program DFG [17] in the complex-pole scheme at next-to- 
next-to-lcading log in QCD contributions and ncxt-to-lcading order (NLO) in electroweak 
contributions. The large parameter space in the MSSM necessitates the use of benchmark 
scenarios [3] . Only the m™o^^ scenario is considered for comparison with previous results. 
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Figure 1: Invariant mass distributions for (a) t^t^, (b) r^Tg, (c) t^t^, (d) t^t/i, (e) TgT^, and (f) 
all candidates. The Z ^ tt background (solid red) is normalised to the theoretical expectation. 
The QCD (horizontal green), electroweak (vertical blue), and Z (solid cyan) backgrounds are 
estimated from data. The ti (vertical orange) and WW (horizontal magenta) backgrounds are 
estimated from simulation and generally not visible. The contribution that would be expected 
from an MSSM signal for M^o = 125 GeV and tan /3 = 60 is shown in solid green. 
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Table 1: Estimated number of events for each background component and their sum, together 
with the observed number of candidates and the expected number of SM signal events for 
Mff = 125 GeV, separated by analysis category. 













TeTh 


Z ^TT 


79.8 ± 5.6 


288.2 ±26.2 


115.8 ±12.7 


146.1 ±9.7 


62.1±8.0 


QCD 


11.7± 3.4 


72.4 ± 2.2 


54.0 ± 3.0 


41.9±0.5 


24.5 ±0.6 


EWK 


0.0 ± 3.5 


40.3 ± 4.3 


0.0 ± 1.3 


10.8±0.5 


9.3 ±0.5 


tt 


<0.1± 0.1 


3.6 ± 0.4 


1.0± 0.1 


< 0.1±0.1 


0.7±0.4 


WW 


<0.1± 0.1 


13.3± 1.2 


1.6± 0.2 


0.2±0.1 


< 0.1±0.1 


Z ^11 


29.8 ± 7.0 






0.4±0.1 


2.0±0.2 


Total 


121.4 ±10.2 


417.9 ±26.7 


172.4 ±13.1 


199.3 ±9.7 


98.7 ±8.0 


Observed 


124 


421 


155 


189 


101 


SM Higgs X 100 


3.9 ± 0.5 


11. 9± 1.6 


3.8 ± 0.5 


9.7±1.3 


4.2 ±0.6 



Both gluon-fusion and associated hh production mechanisms are considered; the former is 
calculated at NLO in QCD using HiGLU [18] with the top-loop corrected to NNLO using 
GGH@NNLO [19], while the latter is calculated at NNLO in QCD using BBH@NNLO [20] 
with the five fiavour scheme. For both SM and MSSM Higgs bosons, the branching 
fractions are calculated using FeynHiggs [21] at the two-loop level. 

The expected distributions of background events are displayed in Fig. 1 and the es- 
timated numbers of events with their associated systematic uncertainties, as well as the 
observed numbers of candidates from data, are given in Table 1. The systematic un- 
certainty on the Z ^ TT background is dominated by the statistical uncertainty on the 
data-driven determination of the efficiency; the other background uncertainties are de- 
scribed in Ref. [12]. 

3 Results 

Limits for model independent and MSSM Higgs boson production are calculated using the 
method of Ref. [22] with CLs = 95% and the test statistic of Eq. 14 from Ref. [23]. The 
test statistic is defined using the profile extended-likelihood ratio of the distributions in 
Fig. 1, where the systematic uncertainties in Table 1 and the uncertainty on the simulated 
invariant mass shapes have been incorporated using normally distributed nuisance param- 
eters. The uncertainty for the invariant mass shape is determined from the momentum 
resolution calibration for simulation, while the primary normalisation uncertainties are 
from luminosity determination and the electron reconstruction efficiency. The distribu- 
tion of this test statistic is assumed to follow the result of Wilks [24] ; this assumption has 
been validated using a simple likelihood ratio. The expected limits have been determined 
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Figure 2: Model independent combined limit on cross- section by branching fraction for a Higgs 
boson decaying to two tau leptons at 95% CLg as a function of M$o is given on the left. The 
background only expected limit (dashed red) and ztlo" (green) and ib2cr (yellow) bands are 
compared with the observed limit (solid black) and the expected SM theory (dotted black) with 
uncertainty (grey). The combined MSSM 95% CLg upper limit on tan/3 as a function of M^o 
is given on the right and compared to ATLAS (dotted maroon and dot-dashed magenta), CMS 
(dot-dot-dashed blue and dot-dot-dot-dashed cyan), and LEP (hatched orange) results. 



using Asimov datasets [23]. 

The upper limit on the cross-section times branching fraction of a model independent 
Higgs boson decaying to two tau leptons with 2.0 < r] < 4.5 is plotted on the left of Fig. 2 
as a function of the Higgs boson mass. The upper-limit on tan /3 for the production 
of neutral MSSM Higgs bosons, as a function of the CP-odd Higgs boson mass, M40, is 
provided in the right plot of Fig. 2. Previously published exclusion limits from ATLAS [25, 
26], CMS [27, 28], and LEP [29] are provided for comparison. 



4 Conclusions 

A model independent search for a Higgs boson decaying to two tau leptons with pseudo- 
rapidities between 2.0 and 4.5 gives an upper bound, at the 95% confidence level, on the 
cross-section times branching fraction of 8.6 pb for a Higgs boson mass of 90 GeV with 
the bound decreasing smoothly to 0.7 pb for a Higgs boson mass of 250 GeV. 

Limits on a MSSM Higgs bosons have been set in the m'j^o^ scenario. Values above 
tan P ranging from 34 to 70 are excluded over the CP-odd MSSM Higgs boson mass 
range of 90 to 140 GeV. For M40 < 110 GeV, these are comparable to the limits obtained 
by ATLAS and CMS using the 2010 data sets but are considerably less stringent than 
the ATLAS and CMS results using 2011 data. The forthcoming running of the LHC 
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should allow the boson, observed by ATLAS and CMS, to be seen in the LHCb detector 
through a combination of channels and should provide complementary information on its 
properties. 
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